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We study the ground state phase diagram of the twisted three-leg spin tube in magnetic fields by the density ma¬ 
trix renormalization group (DMRG) method. The twisted spin tube is composed of triangular unit cells and possesses 
strong quantum fluctuations under geometrical frustration. We apply the sine square deformation method to remove 
strong boundary effects and obtain smooth magnetization curves without steps of finite systems. With the analysis of 
the magnetization curves and correlation functions we determine the ground state phase diagram consisting of (a) a 
Tomonaga-Luttinger (TL) liquid characterized by spin-| Heisenberg model, (b) 3-sublattice state named UUD with 1/3 
magnetization and (c) TL-liquid of massless chirality with 1/3 magnetization plateau, (d) TL-liquid of massless spin 
mode with or without chirality quasi long-range order. 


1. Introduction 

Following Anderson’s proposal of a spin liquid on 2D 
triangular lattice,'* frustrated quantum spin systems have 
been studied to find novel quantum phenomena. Now it is 
widely believed that geometrical frustration leads to degen¬ 
erate ground states and gives rise to diverse low energy prop¬ 
erties of quantum systems. So called ’’chirality” defined by 
cross product of spins on a triangular lattice is a new degree of 
freedom emerged in such systems and the interplay between 
spin and chirality is an interesting research topic of frustrated 
systems. Here we apply external magnetic field and study the 
ground-state properties described by the spin and the chiral¬ 
ity under the magnetic field. Since magnetic field introduces 
unidirectional anisotropy only in spin space, it generally sup¬ 
presses spin fluctuations while leaves the chirality unchanged 
that will lead to a new class of quantum state. 

In this paper we investigate the twisted three-leg spin tube 
consisting of triangular unit cells as a typical ID fully frus¬ 
trated quantum spin system (Fig. 1).^^* The ground state 
of this model has been studied mainly in several limiting 
cases. In the limit of weakly interacting triangles, where intra¬ 
triangle interaction written in the red (thick) lines in Fig. 1 (a) 
are stronger than the other interactions written in the black 
(thin) lines, Fouet et al. have shown that the ground state is 
spin-chirality dimer state with translational symmetry break¬ 
ing.In the opposite limit of strongly interacting triangles, 
they have pointed out that the effective model is a spin-| 
Heisenberg model with a gapless excitation mode. Thus, it 
is expected that the ground state is characterized either by a 
spin-chirality dimer or a spin-| quasi long-range-order, and a 
first order phase transition separates them.^’^* 

In magnetic field, Fouet et al. found 1/3 magnetization 
plateau in the region of weakly interacting triangles.^* When 
the intra- and inter-coupling of the triangles are compara¬ 
ble, Chen et al. have reported that the ground state has 3- 
sublattice structure and UUD state appears in 1/3 magnetiza¬ 
tion plateau. In this plateau, Plat et al. have shown that the 
chirality behaves as XY pseudospin in the limit of weakly in¬ 
teracting triangles"* suggesting the chirality mode is gapless 
while the spin mode is gapful in this limit. Although twisted 
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Fig. 1. (Color online) (a) Structure of twisted three-leg spin tube, (b) 
Equivalent structure of the spin tube with the periodic boundary condition 
in the direction of rung. 


three-leg spin tube is expected to have rich phases, e.g. 3- 
sublattice states, two-component Tomonaga-Luttinger liquid 
of the spin and chirality, and 1/3 magnetization plateau, de¬ 
tailed ground-state phase diagram has not yet been clarified. 

In this paper, we apply the density matrix renormalization 
group (DMRG) method*^* with the sine square deformation 
(SSD)'^* and determine the ground state phase diagram of the 
twisted three-leg spin tube in magnetic field. The paper is or¬ 
ganized as follows. In §2, we define the Hamiltonian of the 
twisted spin tube and explain recently developed technique 
of SSD. In §3 and 4 our numerical results are presented and 
analyzed to determine the ground state phase diagram. We 
summarize our results in §5. 

2. Model and Method 

The Hamiltonian we studied here is defined by 

3 

H = (1) 

1=1 '=1 
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Fig. 2. (Color online) Local magnetization at J\\IJ± = 0-5 obtained by 
DMRG under the usual open boundary condition (OBC) and SSD; (a) lo¬ 
cal magnetization at the center of the system (^^^ 2 ) “ 
site dependence of the local magnetization (5^) — {S\^i 2 ) at = 2.4. Total 
number of unit triangles L is 60 and the typical truncation error in DMRG is 
10“^ 


+ '^11 ^ ^ ■ ^ij+l + ^ij ■ 

f=l ^-=1 

- "lit nr 


>=I '=1 


where Sjj represents spin-1 operator at rung i and leg j, and 
J± (7||) is antiferromagnetic exchange coupling of intra (inter) 
triangles. The last term is Zeeman energy with h being the ex¬ 
ternal magnetic field. To diagonalize this Hamiltonian we use 
the DMRG method, which is usually applied to open bound¬ 
ary conditions (OBC).'^^ However, spins at the ends of the 
open system cause artificial effects called ’’boundary effects” 
that sometimes make it difficult to study the bulk properties 
of the system. Recently, Gendiar et al. introduced the SSD 
and succeed in removing boundary effects. The SSD is a 
kind of energy scale deformation defined as 


L L-\ 

Hssd = X /o0>o(;) + 2 MMiU, j + 1), (2) 

j=i >=i 


where hoij) is on-site term corresponding to the first and third 
term in Eq.(l) and hi{j, j + 1) is the nearest neighbor interac¬ 
tion term such as the second term in Eq.(l). fi(j) is the scaling 
function and defined as 

fi(j) = (3) 

Since the energy scale near the edges of the system is neg- 
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Fig. 3. (Color online) Results at J\\/J± = 2.0 obtained by DMRG: (a) mag¬ 
netization curve (b) local magnetization {S^ ) and (c) entanglement 

entropy S {j, L) of block size j in the L site system under usual OBC. The up¬ 
per (lower) panel in (c) shows the result for = 1/6 (S^/M, = 1/2). 

(d) Dominant local spin structure in the limit of Jy /7x 3> 1.0. 


ligibly small, the boundary effects are efficiently suppressed 
around the center of the system where the scaling function is 
order unity. The bulk properties of the model is then evaluated 
around the center of the system. Eigures 2 (a) and (b) show the 
magnetization curve of {S ^^ 2 ) = 2? (S^ 1 / 2 ) ^i^d the real space 
profile of the local magnetization (5^) = Without 

the SSD, oscillation of (5^) extends into the central part of 
the system and the averaged magnetization is discretized by 
the conservation of the total S^.. In contrast, the oscillation 

J 

under SSD is clearly suppressed and the magnetization curve 
is obtained as a smooth function of h. Such removing of the 
boundary effects and the smooth response to the external field 
are important features of the SSD and these make it easy to 
find anomalies in the response to the external field such as 
magnetization plateaus. The SSD was originally used to re¬ 
store the translational invariance in ID free fermion system 
and it has been shown that the ground state of critical systems 
such as XXZ model and (extend) Hubbard model obtained 
under SSD is identical to the one under the periodic boundary 
condition (PBC).'"'^'®^ In this study we use the SSD for the 
analysis of the ground state correlation functions and magne¬ 
tization curves. 

In our work, we consider the unit triangle as a single site 
and keep up to 1200 basis states in the DMRG calculation 
whose truncation error is the order of 10 ® ~ 10“^. We 
show mainly the numerical results of the finite system of size 
L - 60. Since the numerical results obtained near the edges 
of the system are meaningless in SSD, we only use the results 
around the center. 
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Fig. 4. (Color online) Results at J\\IJ±_ = 1.2 obtained by DMRG: (a) magnetization (^^^ 2 )’ dependence of local magnetization (5? ), (c) transverse 

spin congelation entanglement entropy S{j,L) at = 1/6, 1/3 and 1/2. The upper, middle and lower panels in (b) and (c) 

correspond to /i/Zx = 1.0,1.7 and 3.0, respectively, (e) schematic spin structures of Y, UUD and V states. 


3 . Results 

3.1 7||/yi>1.0 

We start from the case of J\\IJ_i_ > 1.0, where strong inter¬ 
triangle antiferromagnetic couplings ferromagnetically align 
the spins in intra-triangles as shown in Fig.3 (d). The ef¬ 
fective Hamiltonian is then described by spin-| Heisenberg 
model,^' whose low-energy properties are characterized by a 
Tomonaga-Luttinger (TL) liquid. As shown in Fig.3 (a) the 
magnetization at J\\IJ±_ =2.0 monotonically increases with 
the increase in magnetic field and the local magnetization 
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Fig. 5. (Color online) DMRG results at J\\IJ^ = 0.2; (a) magnetization 
curve and (b) entanglement entropy S(j,L) = 1/6 and 1/2. 


{S is site-independent. To confirm this gapless state we cal¬ 
culate the size dependence of the entanglement entropy which 
is analytically obtained for TL liquid as 


S (j, L) = ssln 


L ni 
-sin(-^) 
n L \ 


(4) 


where j and L are the block and the total system size, and sb 
is the constant given by the central charge c as = c/6 (c/3) 
for OBC (PBC). Since the above formula is obtained without 
SSD, we calculate S (j, L) under usual OBC without SSD. Fig¬ 
ure 3 (c) shows S (j, L) at S ^^2 = 1/6 and 1/2, where is 

the saturation magnetization per unit triangle, 3/2. We find the 
central charge c is close to 1 in both cases that indicates the 
ground state is characterized by one-component Tomonaga- 
Luttinger liquid (TLLl) in agreement with the previous work 
by Fouet et al. 


3.2 7||/7x ~ 1.0 

We next consider the intermediate region J\\/J± ~ 1.0. 
Figure 4 (a) shows the magnetization curve at J\\/J± = 1.2 
where we find clear 1/3 magnetization plateau which di¬ 
vides the ground state into three phases; below 1/3 magne¬ 
tization plateau, on the plateau, above the plateau. The up¬ 
per panel in Fig.4 (b) shows the local magnetization 
of the ground state below the plateau. The local magnetiza¬ 
tion {S) j) shows the presence of the 3-sublattice structure. As 
shown in the upper panel of Fig.4 (c), the transverse corre¬ 
lation function (Sf 1 , 2 ^ 71 / 2 + 1 ') 3-sites period. To under- 
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Fig. 6. (Color online) DMRG results in 1/3 plateau at /y/Zx = 0-2; (a) entanglement entropy S(j,L), (b) local magnetization (S^j), (c) local singlet 
correlation (S, j ■ S,+ij) and the z-component of the local chirality ((id)^), (d) longitudinal spin correlation schematic spin structure in 

the 1/3 magnetization plateau at J\\IJi_ < 1.0. In (a), (b) and (d), we also plot the results at J\\IJi. = 1.2 for comparison. The inset in (c) represents local spin 
correlation (S,j ■ Si+ij) for UUD structure at J\\IJ± = 1.2. 


Stand this spin structure let us consider the classical limit. 
When {S ' 11 / 2 +large and positive, xy component of 
the two spins at L/2 and (L/2 + j) sites is parallel, while it 
is almost orthogonal when {Sf^i 2 S~i/ 2 +j'^ ~ antipar¬ 
allel when {S f^ 2 ^ ^ simple picture indicates 

that the ground state below 1/3 plateau has Y structure, which 
is a deformation of 120° state, where two spins are oriented 
upward to the z direction with anitferromagnetic correlations 
in xy components while one spin is antiparallel to z-direction 
as shown in Fig.4 (e). Similarly we can see that the ground 
state above 1/3 plateau is characterized by V structure shown 
in Fig.4 (e). As is shown in the power law decay of the cor¬ 
relation function (‘5,^£/2 ‘^i L/ 2 +/^’ ground states of Y and 
V structures are expected to have massless spin excitations 
consistent with Mermin-Wagner theorem. This is conhrmed 
by the entanglement entropy S (j, L) shown in the upper and 
lower panels of Fig.4 (d), whose j dependence corresponds to 
c = 1 one-component TL liquid. 

In contract to the above results, 1/3 plateau state called 
UUD has short range correlation as shown in the middle panel 
of Fig.4 (c). The local magnetization {S]^ in Fig.4 (b) has 
clear 3-sublattice structure with squeezed moments due to 
quantum fluctuations. The entanglement entropy S (j, L) is in¬ 
dependent of j for large j showing the ground state at 1/3 
plateau is not critical consistent with the short range correla¬ 
tion functions. 

All the 3-sublattice structures Y, UUD and V are origi¬ 
nated from 120° classical state. This was hrst pointed out by 
Chubukov in 2D triangle lattice.They have explained 
these 3-sublattice structures are stabilized by quantum fluc¬ 
tuations. We think the similar 3-sublattice structures in the 


twisted spin tube are also stabilized by the same mechanism 
since the unit cell of the twisted spin tube is identical to that 
of the 2D triangular lattice when J\\IJ±_ - 1.0. 


3.3 7||/7x < 1-0 

We finally investigate the ground state in the region of 
weakly interacting triangles J\\IJi_ < 1.0. Although the mag¬ 
netization curve in Fig.5 (a) shows 1/3 plateau at 7||/7x = 0.2, 
the spin structure in this plateau is different from UUD found 
at J\\IJ±_ ~ 1.0 as will be shown in the following. In the limit 
of weakly interacting triangles, the effective Hamiltonian is 
written as 



+ ^ + h-c.) 



(5) 


where T|(rp represents the raising (lowering) operator act¬ 
ing on the two chirality states |tj = +1/2) and |tj = -1/2) of 
the unit triangle at j-th site, and S j is spin-| operator of the 
same unit triangle.Below the magnetization plateau, 
the bosonization analysis shows that the spin and chirality de¬ 
grees of freedoms are separated within a perturbation analysis 
with respect to the second term in Eq. To conhrm this 

result we first calculate the entanglement entropy S (j, L). As 
seen in the upper panel of Fig. 5 (b), S (j, L) has linear size de¬ 
pendence on In 1^^ sin ^ j and the central charge c is close to 2. 


4 




























































J. Phys. Soc. Jpn. 


FULL PAPERS 


This value of the central charge means the ground state below 
1/3 magnetization plateau is characterized by two-component 
Tomonaga Luttinger liquid (TLL2).^°’^'' Similar result is also 
obtained above the magnetization plateau as shown in the 
lower panel of Fig. 5 (b). We therefore conclude that the el¬ 
ementary excitation above and below 1/3 plateau have two 
massless excitation modes for the spin and chirality. 

We next see /-dependence of S (j, L) in 1/3 plateau. In gen¬ 
eral, the entanglement entropy becomes constant for large / 
if all the excitations from the ground state have a finite gap. 
As seen in Fig.6 (a), S (/, L) at J\\IJi_ - 0.2 has linear depen¬ 
dence on In 1^^ sin with a constant c close to 1. This re¬ 
sult indicates that the elementary excitation has one massless 
mode. Since the spin excitation has a gap in the magnetiza¬ 
tion plateau, the low energy massless excitation of Eq. 5 is 
described by chirality *0 whose Hamiltonian is 

//eft - ^ + h.c.)]. (6) 

,/ 

We therefore conclude that the ground state is characterized 
by TLLl of a chirality XY model consistent with the previous 
work by Plat et a/.'O 

The local spin magnetization (S^j) of the above ground 
state is shown in Fig. 6 (b). In contract to UUD state at 
J\\/J± ~ 1-0, the local magnetization at J\\/J± = 0.2 is uni¬ 
form. The uniform ground state is also shown in the spin cor¬ 
relation function {Sjj-Si+ij} in the same triangle and the local 
chirality expectation value where /F is given by 

3 

= (7) 

1=1 

As shown in Fig. 6 (c), {Sij ■ Si+ij) » -0.25 and ~ 

0.75. These results show the coexistence of a spin singlet and 
chirality in the unit triangle. This is contrasted with the re¬ 
sults at J\\/J± = 1.2 where UUD structure appears with small 
((a^) ) and one positive and two negative local spin correla¬ 
tions, within a triangle. This difference in the 

local correlation makes clear difference in the long range spin 
correlation L/ 2 +j'^ shown in Fig. 6 (d) that is positive 

or nearly 0 at J\\/J± = 0.2, while it has a structure of 3-site 
period with one positive and two negative at J\\/J± = 1.2. 
We illustrate the schematic diagram of the ground state at 
J\\/J± - 0.2 in Fig.6 (e). The ground state is a superposi¬ 
tion of one singlet and one polarized spin like 0.5-0-0 in the 
direction of both rung and leg. Since each polarized spin is 
connected to spin singlet states, there is no penalty in the ex¬ 
change energy among them. Therefor (‘5£/2‘^L/2+j^ 3-site 
period while (5^) is uniform. This is caused by strong geo¬ 
metrical frustration and a unique characteristic of the twisted 
spin tube. 

4. Discussion 

The phase boundaries of the ground state obtained in each 
region J\\/J± > 1.0, J\\/J± 1.0 and J\\/J± <1.0 are pre¬ 

sented in the phase diagram shown in Fig.7 (a). The boundary 
is determined by the appearance or the change of 3 sublat¬ 
tice structure and the central charge of TL liquids. For the 
boundary between V and TL liquid states in the region of 
J\\/J± > 1-0, we have used the appearance of the 3 sublattice 
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Fig. 7. (Color online) (a) Ground state phase diagram of the twisted spin 
tube. The labels D, TLLl, TLL2 represent the dimer state, one- and two- 
component Tomonaga-Luttinger liquid, respectively. The labels Y and V rep¬ 
resent Y, V states. PuuD and Pch corresponds to 1/3 magnetization plateau 
with UUD structure and 1/3 plateau with TLLl of chirality, respectively. FP 
is fully polarized state, (b) local magnetization at J\\/J± = 1.4 as a 

function of /i/Zx- (c) entanglement entropy S {j, L) = 1/3. 


Structure as a signal of the transition. This is clearly shown in 
the local spin polarization - 1-4- As shown 

in Fig.7 (b), splits into two values at /i r; 4 which 

corresponds to the change from spin-| TLLl to V state. The 
phase diagram around the plateau state is rather complicated. 
Since the corresponding 2D triangular lattice does not have 
1/3 plateau in the region of J\\IJ± < 0.5,the appearance of 
the magnetization plateau with massless mode of chirality is a 
unique feature of the twisted spin tube. Above and below the 
magnetic plateau, the elementary excitations at 7||/7± < 0.3 
are characterized by the central charge c ^2 and are different 
from 3-sublattice state and TLLl. 

Ligure 7 (c) shows S(j,L) for different J\\/J± in 1/3 mag¬ 
netization plateau. It is seen that S (j, L) has linear size de¬ 
pendence on In sin at 7i|/7j_ = 0.18 and 0.20, while it 
dose not at - 0.26. Thus, TLLl of the chirality with 

uniform magnetization extends to 7||/7j_ ~ 0.26. We expect 
the transition from TLLl of the chirality to UUD state is sec¬ 
ond order because we find no signal of level crossing, but it is 
difficult to confirm this point within a finite system. Detailed 
finite size scaling analysis is needed to clarify the nature of 
the transition. 

5. Summary 

We have studied the ground state of the twisted three-leg 
spin tube in magnetic field by the DMRG with the SSD and 
identified various phases such as TLLl, TLL2, 3-sublattice 
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state and 1/3 magnetization plateau. In particular we have 
confirmed that for J\\IJi. < 0.3 the ground state in the 1/3 
plateau has low energy excitations described by massless 
mode of chirality. The coexistence of the spin and chirality 
degrees of freedom and the interplay between them are char¬ 
acteristic features of twisted three-leg spin tube and the origin 
of its diverse ground states. 
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